Multiferroic behavior of Aurivillius Bi 4 Mn 3 Oi2 from first-principles 
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The multiferroic behavior of the hypothetical Aurivillius compound Bi4MnaOi2 has been explored 
on the basis of density functional calculations. We find that the tetragonal paraelectric phase of this 
material is ferromagnetic, showing ferroelectric and antiferrodistortive instabilities similar to the 
ones observed in its ferroelectric parent compound B14T13O12. Our results indicate, however, that 
the presence of Mn +4 ions at the B-sites shrinks the cell volume and consequently the unstable polar 
mode, associated with the ferroelectric polarization, is overcame by an antiferrodistortive distortion. 
In this way, Bi4Mn30i2 exhibits incipient ferroelectricity at its equilibrium volume. We show that 
the ferroelectric state can be favored by strain or partial substitution of Mn with Ti. 

PACS numbers: 77.84.-s, 75.85.+t, 63.22.Np 
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I. INTRODUCTION 

Multiferroic materials, exhibiting both ferroelectric- 
ity and magnetic order in the same phase, are of par- 
ticular interest for their potential technological appli- 
cations. One of the current trends in the search for 
new multiferroic compounds is, starting from a ferroelec- 
tric host, to incorporate magnetically active species and 
check whether it is (anti)ferromagnetic and insulating. 
An alternative way to predict multiferroic behavior is 
from material computational design, that is the predic- 
tion of multiferroic properties using first-principles tech- 
niques prior to their experimental verification; 1 .^ Exam- 
ples of this methodology include the prediction of a ferro- 
electric instability in BiMnOsj 4 - several unstable phonon 
branches in BiCr03 p structural instabilities of CaMnOs 6 
and the properties of as-yet-unsynthesized 3d-5d double 
perovskites.— 

Between the ferroelectric materials, Aurivillius layered 
perovskites are good trials for multiferroics. Described 
by the formula [Bi202][A n _iB n 03„+i], they are formed 
by stacking fluorite blocks of Bi202 with n perovskite- 
like blocks. Various A and B cations are allowed, so 
that this family has numerous representatives. Archetype 
compounds are Bi2WOg (n=l), SrBi2Ta20g (n=2) and 
Bi 4 Ti30i2 (n=3). Recently some experimental groups 
have synthesized Aurivillius compounds incorporating 
magnetic ions at the i?-site. The underlying idea was 
the substitution of magnetic transition metal cations 
into the central octahedron layer of the perovskite-type 
block, while maintaining ferroelectric displacements in 
the outer octahedron layers. For that reason, those works 
involved high-order Aurivillius phases (n>3), such as 
Bi 5 Ti 3 FeOi 5 ^ Bi 6 Ti 3 Fe 2 18 ^ and Bis^CrOig. 11 
Sharma et al. reported ferroelectric behavior in a fam- 
ily of transition metal substituted three-layer Aurivil- 
lius phases, formulated Bi2_ K Sr2+a(Nb/Ta)2+ a ;Mi_ a :Oi2 
(x w 0.5, M = Ru 4 +, Ir 4 + or Mn 4 +), showing that 



these compounds have the same orthorhombic symmetry 
and similar dielectric and ferroelectric properties as their 
(non-magnetic) ferroelectric parent compounds j 12 : 13 Sev- 
eral investigations have been based on the incorpora- 
tion of perovskite like LaFeOs, BiMn03, and BiFe03 
into a Bi 4 Ti30i2 matrix. 14-16 However, one can find that 
the reported magnetic properties in these materials are 
quite weak even at low temperature. The modification of 
B14T13O12 ceramics was also investigated by substituting 
their cations partly with Gdjil Fe^ and Mn 4^ 

In spite of the large amount of experimental work 
showing that the generic structure of Aurivillius phases 
has the potential to serve as a template for multiferroic 
behavior, theoretical studies that help to clarify the ef- 
fect of magnetic species in the Aurivillius structure are 
missing. In this paper we investigate the effects of the 
substitution of Ti ions by magnetically active Mn +4 ions 
at the B sites of Bi4Ti30i2 (BIT), which is a proto- 
typical Aurivillius compound with bismuth also occu- 
pying the perovskite A site. BIT has received tremen- 
dous attention during the last decade, especially in the 
thin film form, due to its promise in nonvolatile ferroelec- 
tric random access memory applications.— It has a rea- 
sonably large spontaneous polarization (~ 50 ^C/cm 2 ), 
a high Curie temperature (950 K), and fatigue resis- 
tance on Pt electrodes. Its complex crystal structure can 
be described in terms of relatively small perturbations 
from a high-symmetry body-centered tetragonal struc- 
ture (space group symmetry I4/mmm), which contains 
only one formula unit per primitive cell. Two main dis- 
tortions from this tetragonal (paraelectric) structure lead 
to the ferroelectric orthorhombic phase. First, the ions 
displace along the [110] axis of the tetragonal structure. 
Second, the TiOe octahedral rotate around the a and c 
axesi^i The first distortion is directly responsible for the 
observed macroscopic polarization along the orthorhom- 
bic a direction ([110] axis of the tetragonal structure). 

Specifically, we performed density-functional theory 
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calculations to investigate the hypothetical three-layer 
Aurivillius compound Bi4Mn 3 0i2 (BIM), which results 
from substituting the three _B-site Ti 4+ cations with 
Mn 4+ in the BIT lattice. We predict structural, elec- 
tronic and magnetic properties of BIM, deepening the 
understanding of the effects of magnetic species in the 
Aurivillius structure of BIT. The paper is organized as 
follows. In Sec. II we describe computational details. 
Results are presented in Sec. III. Within this section we 
compare our BIM results with those of BIT from previ- 
ous works and also with results we have obtained. First, 
the tetragonal structure of BIM is optimized and the elec- 
tronic structure and magnetic order of this phase are dis- 
cussed. Then we focus on the two more unstable modes 
of BIT: the zone-center E u modes related to the polar 
instability and the zone-boundary modes related to 
oxygen antiferrodistortive (AFD) distortions. We inspect 
then the energy landscape along the unstable modes di- 
rections and hnd that by replacing the Ti cations with 
Mn, the polar E u mode is overcome by the Xg~ mode. 
Finally, we investigate the possibility of favoring the fer- 
roelectric distortion by strain and chemical engineering. 



II. COMPUTATIONAL DETAILS 

We performed standard first-principles calculations us- 
ing density-functional theory within the generalized gra- 
dient approximation (GGA) with the Perdew, Burke 
and Ernzcrhof 2 '' parametrization as implemented in the 
VASP-4.6 code^i The projector augmented wave poten- 
tials were used including explicitly 15 valence electrons 
for Bi (5d 10 6s 2 6p 3 ), 13 for Mn {3p 6 , 3d 5 , 4s 2 ), 10 for 
Ti {2p 6 , 3d 3 , 4s 1 ) and 6 for oxygen (2s 2 , 2p 4 ). Con- 
vergence with respect to the fc-point mesh was achieved 
using a 6x6x6 T-centered grid. For a proper treatment 
of strongly localized Mn d electrons, we used the on-site 
Coulomb correction (GGA+J7) in the Dudarev imple- 
mentation with parameter U e ff = U — J = 5.2 eV£l 
Except for the calculations done to determine the mag- 
netic ground-state, where different U e // values were used 
to check the energy differences between spin configura- 
tions, all other calculations were performed using the 
mentioned value. 

Using a 19-atom body-centered-tetragonal unit cell 
with space group I4/mmm, we relax all degrees of free- 
dom (lattice parameters and internal coordinates along 
the stacking direction z) for spin polarized structures. 
We optimized the ground state structures with a plane- 
wave cutoff of 600 eV which yielded good convergence 
and, hence, was also used for the rest of the calcula- 
tions. To search for the presence of ferroelectric and 
antiferrodistortive instabilities in the tetragonal para- 
electric structure, we determined the phonon frequencies 
and eigenvectors of the infrared-active E u and rotational 
X$ modes, which are the main unstable modes of BIT. 
To this end, we calculated atomic forces of symmetry- 
adapted configurations determined with the help of the 




FIG. 1. (Color online) BIM tetragonal phase along [001] from 
to c/2. 



ISOTROPY software^ where selected atoms are dis- 
placed by 0.01 A along the x, y and/ or z direction. To 
compute the X$ modes, in particular, the primitive cell 
was doubled to a 38-atom conventional tetragonal cell. 
From the forces as a function of displacements, a dy- 
namical matrix of 11x11 (8x8) for E u (X3 ) symmetry 
was constructed and diagonalized. Once the eigenvec- 
tors have been determined, we evaluate the total energy 
as a function of the displacement pattern of the unstable 
(imaginary frequency) modes. Those curves will provide 
the ferroelectric and antiferrodistortive instability energy 
associated with a particular phonon mode. 



III. RESULTS 

We begin by describing the optimized tetragonal struc- 
ture of BIM. The minimized lattice parameters are a = 
3.811 A and c = 32.613 A. These values are smaller than 
those we compute for BIT using GGA: a = 3.842 A and 
c = 33.241 A, which in turn differ slightly from the 
experimental values reported in the literature (a = 
3.8633 A and c = 33.2942 A^ and a = 3.850 A and 
c = 32.832 A)2i. As seen the replacement of Ti 4+ ions 
with Mn 4+ shrinks the equilibrium volume of the unit 
cell, which is consistent with the fact that six coordi- 
nated Mn 4+ has an ionic radius (0.67 A) smaller than 
that of Ti 4+ (0.745 A). This chemical substitution would 
be analogous to applying pressure on BIT and one could 
expect, following the general trend in many perovskite 
oxides, a strong sensitivity of the ferroelectric instability 
to volume. 

In the perovskitelike blocks of the tetragonal structure, 
by symmetry there is one central Mn 4+ ion [Mn(l) in 
Fig.[Tj and two equivalent Mn 4+ ions above and below the 
central layer, Mn(2). In the optimized structure of BIM, 



TABLE I. Wyckoff positions for relaxed BIM tetragonal struc 
ture are compared to GGA calculations and experimental re 
suits for BIT. 



Atom 


Site 


BIM 


BIT 


BIT cxp. 29 


Bil 
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0.4342 


0.4336 


0.4325 


Bi2 
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0.2901 


0.2887 


0.2885 


Mnl 
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0.0000 


0.0000 


0.0000 


Mn2 
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-0.1254 


-0.1275 


-0.1293 


Ol 
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0.00000 


0.0000 


0.0000 
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0.3858 


0.3848 


0.3823 


03 
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-0.0593 


-0.0595 


-0.0589 


04 
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-0.1810 


-0.1815 


-0.1817 


05 


d 


0.0000 


0.0000 


0.0000 



the two Mn(2) cations are pushed away along the c-axi: 
from the center of each oxygen octahedron towards thi 
fluorite blocks. Besides, the oxygen octahedral arount 
these Mn(2) ions are considerably elongated along the c 
axis, opposite to the undistorted central octahedron. Th< 
A-site Bi 3+ ions [Bi(l)] also move towards the fiuoriti 
blocks, that is, they approach the 0(2) ions, moving awa^ 
from the 0(1) ions. Table |TJ lists the optimized interna 
atomic positions of BIM. For comparison, we include ou: 
GGA results for BIT and the respective experimenta 
values from Ref. H^. As seen the internal coordinates ii 
both compounds are similar. 

The magnetic order of BIM was investigated by in- 
specting the total energy of different configurations with 
spin pointing along [001]. We have optimized separately 
the tetragonal structure for the ferromagnetic (FM) or- 
dering and for the following antiferromagnetic (AFM) 
arrangements in the perovskitelike block: A, C, and 
G-type, in which the fraction of parallel aligned first- 
neighbor spins is 2/3, 1/3, and 0, respectively. By com- 
paring the total energy of those configurations we find 
that the magnetic ground state of the material is FM 
with local magnetic moments of 3.3 and 3.4 /j,b for Mn(l) 
and Mn(2), respectively, which are slightly higher than 
the formal value of 3 \ib expected for a four- valence Mn 
(d 3 ) atom. We find that the A- type AFM arrangement is 
located at an energy slightly higher than the FM state, 
just a few mRy/f.u. above it. The energy of the C-type 
AFM ordering is more separated, ~ 15 mRy/f.u. above 
the FM state. Finally, the G-type AFM configuration 
lies a few mRy/f.u. above the G-type. It is known that 
the energy differences between spin configurations are di- 
rectly related to the magnetic interactions or exchange 
constants, which depend on the on-site Coulomb repul- 
sion parameter U e //. Then, we have verified that the FM 
order is the most stable configuration for different choices 
of U e / f . Although the variation of U e / / produced quan- 
titative changes, the hierarchy of the magnetic configu- 
rations is kept. We point out that the optimized crystal 
structure of BIM described above corresponds to the FM 
ground state. Furthermore, all of the following results 
have also been obtained using the FM spin configuration. 

The total density of states (DOS) and projected DOS 
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FIG. 2. Total majority-spin state and projected DOS onto 
atomic sites for tetragonal BIM. Zero energy is set as the 
Fermi level. 



onto the atomic sites for the majority spins are shown 
in Fig. [2] The system has no band gap in the paraelec- 
tric I4/mmm phase, similar to what was observed for 
the paraelectric phase of other ^4Mn03 perovskites.— ^ 
Instead, a half-metallic electronic structure is achieved 
with a band gap of ~ 2.2 eV in the minority spins (not 
shown here) . The occupied states of the conduction band 
(ranging from -7 cV to the Fermi level) are mainly O 2p 
and Mn 3c? states, which are strongly hybridized. The 
low-energy Bi states, around -10 eV, come from the 6s 
orbitals (the so-called lone pairs), while the unoccupied 
Bi 6p states are located at high energies around 6 eV. 
The contribution of the Bi 6s and 6p states to the con- 
duction band indicates that these states are hybridized 
with the 2p of their neighbor oxygens, which is similar to 
what has been reported for BIT. 22 To better visualize the 
differences between BIM and BIT, the partial DOS of the 
-B-site ion which is in an octahedron crystal field, are pro- 
jected onto the t2 g and e g orbitals (black lines and gray 
dashed areas in Fig. G2 respectively). The DOS of BIT 
was obtained at the BIM equilibrium crystal structure to 
avoid differences that could arise from the different lat- 
tice parameters of the two compounds. It is clear from 
the figure that the majority occupied Mn 3d states are 
much more extended than those of Ti, showing that the 
Mn 3d hybridize more strongly with the O 2p states. The 
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FIG. 3. DOS projected onto £29 (black curves) and e 9 (gray 
dashed regions) orbitals for each Mn site in BIM for the ma- 
jority and minority spins. They are compared to the respec- 
tive orbitals of the Ti cations in BIT. The DOS of BIT was 
computed at the BIM equilibrium volume. 



majority-spin Mn ti g states are fully occupied whereas 
the corresponding minority spins are practically empty. 
Concerning the e g orbitals, the majority spins are more 
extended and they cross the Fermi level, which is respon- 
sible for the paraelectric phase being half-metallic. 

Next we investigate the possibility of ferroelectricity 
in BIM by studying the unstable normal modes of the 
high-symmetry I4/mmm tetragonal phase, and compare 
our results with those for BIT. We note that as BIM 
is a hypothetical compound, we do not have an exper- 
imental low-temperature ground-state structure to de- 
compose and identify the symmetry of active modes that 
would freeze in the high-temperature tetragonal phase. 
That limitation leads us to analyze the normal modes in 
tetragonal BIM with the same symmetry of the relevant 
modes found in the tetragonal parent compound BIT. 
Furthermore we are aware that the energy landscape of 
Aurivillius compounds, as it was shown in BIT 23 and 
SrBi2Ta20g, 24 is much more complex than that of typical 
ferroelectric oxides, in which the freezing of one unsta- 
ble polar mode directly explains the ferroelectric phase. 
So, based on an exhaustive study of the BIT phonon 
modes by Perez-Mato et. alj£. we restrict our discussion 
to zone center polar modes with E u symmetry and zone 
boundary antiferrodistortive modes with symmetry. 
We find that two of the 11 bi-dimensional E u modes of 
BIM have imaginary frequency, similarly to what was ob- 
tained for BIT The most unstable E u mode [E u (l), 
uj = « 146 cm -1 ] involves primarily antiphase displace- 
ments of the Bi(l) ions with respect to the apical oxygen 
0(3), plus a smaller displacement of the central Mn(l) 
against its neighbor oxygens 0(1). A similar pattern 
was indeed associated with the in-plane ferroelectric po- 



larization of BITi 2 ^ The second unstable mode [E u (2), 
uj = i57 cm ] can be roughly described by the displace- 
ments of the Bi202 layers relative to the perovskitelike 
blocks. This kind of instability was found in other Au- 
rivillius family memberSj 31 and seems to be a general 
feature of them. Surprisingly, not only the eigenvectors 
but also the frequencies of the two unstable E u modes 
are very similar to the ones obtained in a linearized aug- 
mented plane-wave (LAPW) calculation of BIT, where 
uj = il45 cm and uj = i48 cm -1 were assigned for 
the E u (l) and E u {2) modes respectively^ Regarding the 
zone-boundary modes, the most unstable mode [Jfg (1), 
uj = «242 cm -1 ], is much more negative than the po- 
lar ^u(l) mode. This mode takes place mainly at the 
center of perovskite slabs, tilting the central octahedron 
around the [110] direction, a~a~c° in the Glazer nota- 
tion. However, as the outer apical 0(4) oxygens almost 
do not move, the two external octahedral are distorted. 
The second unstable AFD mode [X % (2)] consists mostly 
of the tilting of the two external octahedral around the 
[110] direction. This mode involves the displacement of 
the 0(4) oxygens as well as the Bi(2) ions belonging to 
the fluorite layer. In summary, the described eigenvectors 
of the four unstable E u and X J modes are very similar 
to those obtained in BIT (illustrations of the different 
unstable eigenvectors of BIT are included in Ref. 23i). 

The energy landscape around the tetragonal structure 
as a function of the amplitude in A of the unstable modes 
is shown in Fig. Hta). For each mode, the eigendisplace- 
ment vector is normalized with respect to the tetragonal 
high-symmetry 38-atom unit cell. We summarize our re- 
sults as follows. First, the most unstable A^"(l) mode 
decreases the total energy in ~ 21.4 mRy/f.u., against 
the ~ 7.5 mRy/f.u. produced by the E u (l) polar mode. 
This is strikingly different from what is observed in BIT 
where the rotational mode is overcame by the in-plane 
polar instability. With regard to the other two modes, 
X% (2) and E u (2), they present a rather small energy 
instability. However, when each one by separate is com- 
bined with the most unstable mode of the same symmetry 
(open symbols), they further decrease the energy mini- 
mum of the double- well obtained with the most unstable 
mode. Interestingly, the particular an-harmonic coupling 
of X^ (2) with X^(l) (open squares) yields an absolute 
energy minimum. The combination of these two rota- 
tional modes results in an atomic structure that differs 
from the high-symmetry phase mainly in that the cen- 
tral octahedron is slightly distorted and rotated around 
the [110] direction at an angle of ~ 10°. This structure, 
however, is still half-metallic and ferromagnetic, as is the 
high-symmetry one. Finally, in Fig. 4(b) and 4(c), we 
explore the an-harmonic coupling between the two most 
unstable modes of different symmetry, X3 (1) and E u (l), 
to investigate whether they cooperate in the energy min- 
imization. We find that the AFD Xg~(l) distortion fully 
stabilizes the E U (X) polar mode [open circles in Fig. 4(b)], 
when the last one is computed from the structure corre- 
sponding to the energy minimum of X^ (1), compared to 
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FIG. 4. (Color online) (a) Total energy per formula unit rela- 
tive to the optimized tetragonal structure as a function of the 
unstable mode amplitudes with symmetry: (filled trian- 
gles and squares) and E u (filled circles and diamonds). Com- 
bination of same symmetry modes: less unstable mode X^(2) 
(open squares) and E u (2) mode (open diamonds) are super- 
imposed to the energy minima of the most unstable mode, 
X^(l) and (E u (l)), respectively, (b), (c) Coupling effect be- 
tween the two most unstable modes of different symmetry, (b) 
Change of energy per f.u. as a function of the E u (l) mode am- 
plitude relative to a structure with mode X3 (1) frozen with 
amplitude of 1 A (open circles). The curve for the pure E u (l) 
distortion (filled circles) is repeated for comparison with, (c) 
Change of energy per f.u. as a function of the Xf(l) mode 
amplitude relative to a structure with mode E u (l) frozen with 
amplitude of 0.8 A(open triangles), and for the pure Xg"(l) 
distortion (filled triangles). 



TABLE II. BIM relaxed lattice parameter in the tetragonal 
structure at two different in-plane strains and in the equilib- 
rium value. 



In-plane strain 


a (A) 


c(A) 


-2 % 


3.7345 


33.5498 





3.8108 


32.6128 


+2 % 


3.8870 


32.0985 



the pure E u (l) distortion from the optimized tetragonal 
structure. Instead, the polar mode effect on the AFD 
X^(l) [open triangles in panel (c)] is not so strong, but 
it does reduce the AFD double-well instability. There- 
fore, these behaviors indicate that the coupling of 
and £^(1) is penalized. 

Based on the existence of the polar instability in tetrag- 
onal BIM (although it is overcame by the rotational dis- 
tortion), we investigate to what extent it would be pos- 
sible to favor ferroelectricity by the application of epi- 
taxial strain. We analyze two different structures around 
the equilibrium volume with in-plane lattice parameter, 
a + 2% and a_2%- The lattice parameters obtained by fix- 
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FIG. 5. (Color online) Effect of epitaxial strain on the two 
most unstable normal modes: E u (left panel) and X^ (right 
panel). Energy differences relative to the optimized tetrago- 
nal structure for a eq (filled lines), a +2 % (dashed lines), and 
a_2% (dotted lines) using the eigenvectors obtained at the 
equilibrium volume. Open symbols correspond to BIT at the 
BIM equilibrium volume. 



ing the in-plane lattice constants, and fully relaxing c and 
the atomic positions, are listed in Table [Til 

Figure [5] displays the energy landscape of the strained 
structures and compares them with the stress-free one. 
As expected, the in-plane E u mode (left panel) is favored 
by tensile bi-axial strain, whereas the AFD instabil- 
ity (right panel) is disfavored. This behavior, which is 
typical in ABO3 perovskite oxides, indicates that epi- 
taxial strain seems to be an appropriate route to induce 
ferroelectricity in this Aurivillius compound. However, it 
can be estimated by extrapolating the above data that 
it is necessary to apply a strain of ~ +6% to favor the 
polar mode over the antiferrodistortive distortion. This 
crudely estimated value is out of the range of strains usu- 
ally accessible through epitaxial growth. 

Since applying a large epitaxial strain is not easily 
achievable in practice, we instead explore the possibility 
of making BIM ferroelectric under a partial substitution 
of Mn by Ti. As a first test, we take the BIM equilib- 
rium structure and replace all of the Mn ions with Ti; 
that is, we analyze BIT at the BIM equilibrium volume. 
The resulting energy curves are also plotted in Fig. [3] 
with open symbols. By comparing them with the BIM 
curves (filled lines), we see that effectively the presence of 
the Ti ions favors the polar instability and considerably 
decreases the AFD one. We can then speculate that a 
partial substitution of Mn with Ti would also favor fer- 
roelectricity. To support this idea we check the energetics 
of the Aurivillius compounds Bi4(TiMn2)Oi2 with the B 
sites in the perovskitelike layer occupied by one Ti at 
the central octahedron, and two Mn ions in the external 
ones. We set the lattice parameters a and c by interpo- 
lating those values between the optimized BIM and BIT 
ones. We find that the polar mode E u indeed is more 
unstable than the AFD one by ~ 3.8 mRy/f.u., which 
shows that the Ti ion at the central octahedron plays an 
important role in stabilizing the ferroelectric mode over 
the AFD one. This is consistent with the fact that the 
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two modes of interest, as described above, are localized 
in the center of the perovskite slab. We think that this 
particular solid solution deserves a more detailed study 
in the future. Nevertheless, we can conclude that our 
results indicate that, not only the cell volume, but also 
the nature of the transition metal ion is indeed relevant 
for the stabilization of the ferroelectric phase in Aurivil- 
lius compounds. Therefore, it would be also interesting 
to further investigate the effect of magnetic species with 
higher ionic radius than Mn 4+ , such as Ru 4+ and V 4+ . 

In summary, we have investigated from first-principles 
the effects of magnetic Mn species in the Aurivillius struc- 
ture of Bi4Ti30i2- Specifically we explored the possibil- 
ity of finding multifcrroic behavior in the hypothetical 
Aurivillius compound Bi4Mn 3 0i2, which results from re- 
placing the Ti with Mn ions. We have found that the 
paraelectric tetragonal structure of this material presents 
ferromagnetic order and structural instabilities similar to 



BIT. We have shown, however, that the presence of Mn 
ions shrinks the cell volume, and the polar mode associ- 
ated with the in-plane ferroelectric polarization, is over- 
came by an antiferrodistortive distortion. In this way, 
Bi 4 Mn 3 0i2 exhibits incipient ferroelectricity at its equi- 
librium volume. We have seen that although it is possi- 
ble to favor the polar state by the application of epitaxial 
strain, the necessary strain value is not easily achievable 
in practice. Our results indicate that the partial sub- 
stitution of Mn with Ti is a fruitful alternative to favor 
the ferroelectric distortion, highlighting the relevant role 
played by the nature of the transition metal ion. We thus 
hope that our findings may be useful to better understand 
the multiferroic properties of Aurivillius compounds. 
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